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Abstract 

The Ca 3 and Ca 4 metallic clusters have been investigated using state-of-the art ab 
initio quantum mechanical methods. Large atomic natural orbital basis sets have been used 
in conjunction with the singles and doubles coupled-cluster (CCSD) method, a coupled- 
cluster method that includes a perturbational estimate of connected triple excitations, 
denoted CCSD(T), and the multireference configuration interaction (MRCI) method. The 
equilibrium geometries, binding energies and harmonic vibrational frequencies have been 
determined with each of the methods so that the accuracy of the coupled-cluster methods 
may be assessed. Since the CCSD(T) method reproduces the MRCI results very well, cubic 
and quartic force fields of Ca 3 and Ca 4 have been determined using this approach and used 
to evaluate the fundamental vibrational frequencies. The infrared intensities of both the 
e' mode of Ca 3 and the t 2 mode of Ca 4 are found to be small. The results obtained in this 
study are compared and contrasted with those from our earlier studies on small Be and 

Mg clusters. 
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Introduction 


There has been considerable recent interest in the properties of small clusters (see 
for example Refs. 1-10), motivated principally by two issues. The first is the question 
of convergence of cluster properties towards the bulk values. Of course some properties 
will approach the bulk value more quickly than others as the cluster size is increased. The 
second issue is interaction between theory and experiment. The study of small clusters has 
progressed very rapidly since accurate experimental studies may be used to evaluate the 
predictive reliability of different theoretical methods, and then accurate theoretical studies 
may be used to evaluate or aid in the design of new experimental techniques. 

Our studies of small clusters 11 ’ 12 have focused on computing the structures, binding 
energies, vibrational frequencies and infrared intensities of the trimers and tetramers of the 
alkaline-earth elements beryllium and magnesium using elaborate treatments of electron 
correlation. We have also examined 13 the equilibrium structures and binding energies of 
the pentamers, Be s and Mg 5 . In the present study we extend our investigations to in- 
clude Ca 3 and Ca 4 . There are three previous studies 14 ’ 15 of Ca 4 (and none of Ca 3 ) that 
have incorporated electron correlation effects. In one Ca 4 study a total binding energy of 
18.3 kcal/mol, relative to four Ca atoms, was obtained at the single-reference single and 
double excitation configuration interaction (CISD) level of theory, including Davidson s 
correction 16 for higher excitations. In the other, a binding energy of xx.x kcal/mol was ob- 
tained using a multireference Cl approach (based on SCF orbitals). Based on our previous 
studies of the Be and Mg clusters, it is likely that even the higher value is a substantial 
underestimate of the true binding energy. 

Another important component of our earlier studies 11-13 ’ 17 is the comparison of ge- 
ometries and binding energies obtained with various electron correlation methods. The 
s-p near-degeneracy effects in the alkaline-earth valence shell are very large, and strongly 
influence the binding in the small clusters. It is not only essential to describe these non- 
dynamica! effects accurately, however, but also to account properly for dynamical correla- 
tion in order to obtain reliable binding energies for these systems. Hence the most desirable 
treatment might appear to be a full valence complete-active-space SCF (CASSCF) calcu- 
lation followed by a multireference configuration-interaction (MRCI) calculation. This is 
indeed an excellent level of treatment, but unfortunately it becomes very expensive to 
apply in large basis sets and at many geometries. 
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In our studies of the lighter alkaline- earth clusters, we have made extensive use 
of the single and double excitation coupled-cluster approach (CCSD) corrected with a 
perturbational estimate of triple excitations (CCSD(T)). 18 The CCSD(T) method performs 
very well in comparison with MRCI results for the lighter alkaline-earth clusters, and 
appears to treat both the non- dynamical and dynamical correlation effects in these systems 
accurately. We have also established the reliability of the CCSD(T) method by full Cl 
comparisons on Be 3 . 20 compared results obtained at the MRCI level of treatment with 
those obtained using the single and double excitation coupled-cluster (CCSD) approach 
and the CCSD(T) method, in which a perturbational estimate of the effects of connected 
triple excitations has been included. 18 In these studies we established that the CCSD(T) 
method performs very well in reproducing results obtained from the MRCI approach. 
Based on many comparisons of full configuration interaction (FCI) and MRCI energies, 1 
it has been asserted that the MRCI results should be very close to the n-particle limit, 
and we have established 20 that the MRCI and CCSD(T) results for Be 3 are indeed very 
close to the analogous FCI quantities. For comparison purposes, we have again used the 
CCSD, CCSD(T) and MRCI methods in examining the potential energy surfaces of Ca 3 
and Ca 4 so that it may be determined how well the coupled-cluster methods perform for 
the series of alkaline-earth metals. 

In the next section we describe the computational methods employed in this study 
and in the following section our results are presented and discussed. A comparison of our 
new results for the Ca clusters and our previous results for the Be and Mg clusters is also 
given. The final section contains our conclusions. 


Computational methods 

Two atomic natural orbital 21 (ANO) basis sets have been used in this study. The 
(22s 17p) primitive basis set is that of Partridge 22 and was augmented with a (Ad 3/) even 
tempered polarization set defined by a = 2.5 n o:o for n = 0,...,k. The 0=0 values for the d 
and f functions are 0.0232 and 0.0440, respectively. The smallest basis set consists of 5s, 

4p and Id ANOs and will be denoted [5s 4 p Id}. The larger basis consists of 6s, 5 p, 2d 
and If ANOs and will be denoted [6s 5p 2d If}. The ANO contraction coefficients were 
obtained by averaging the natural orbitals from single and double configuration interaction 
(CISD) calculations on the lowest and l P states of atomic Ca. Only the pure spherical 
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harmonic components of the d and / functions have been used. 

As discussed in the Introduction, the CCSD, CCSD(T) and MRCI methods have 
been used to treat electron correlation. In all cases, only the Ca 4s electrons have been 
included in the correlation procedure. The coupled-cluster wave functions are based on 
self-consistent field (SCF) molecular orbitals while the MRCI wave functions are based 
on CASSCF molecular orbitals. All valence electrons were allowed variable occupancy m 
all valence orbitals in the CASSCF calculations. References for the MRCI wave functions 
were selected using a 0.05 threshold — that is, all occupations having a component spin- 
coupling with a coefficient of 0.05 or larger in the CASSCF wave function were used as 
references in the MRCI procedure. 

In analogy with small Be and Mg clusters, the equilibrium geometries of Ca 3 and 
Ca 4 were constrained to have D 3h and T d symmetry, respectively. Harmonic frequency 
analyses demonstrate explicitly that these geometries are indeed minima on the Ca 3 and 
Ca 4 potential energy surfaces. In addition, a linear structure for Ca 3 was optimized and 
found to be significantly higher in energy than the equilateral triangle. Hence it is expected 
that the equilateral triangle and tetrahedron are the global minima on the Ca 3 and Ca 4 
potential energy surfaces, respectively. 

The quadratic, cubic and quartic force constants of Ca 3 and Ca 4 have been deter- 
mined numerically and are given in symmetry internal coordinates. The symmetry internal 
coordinate definitions are: 


Ca 3 


Ca 4 


Si (4) = + r 2 + r 3) 

S 2 a(e') = ^(2r!-r 2 -r 3 ) 

s 2b(e) = -j={r 2 -r 3 ) 


5i(a a ) 

S 2 a(e) 

S 2 b{e) 


^(r 1 +r a +r,+r 4 + r,+r 6< 


— ^=(2r 3 - r 2 - r 3 + 2r 4 — r 5 - r 6 ) 
Vl2 


— — (r 2 — r 3 + r 5 — r 6 ) 

V4 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 
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s 3 x {h) = ^ 

Szy{h) = ^( 7 ’ 3 " 7 ’ 6 ) ^ 

5 3 z(i 2 ) - 7f( ri_r4 ^ ^ 

For Ca 4 , the numbering of bonds is such that if n connects one pair of atoms, r 4 connects 
the other pair, and similarly for the pairs {r 2 ,r 5 } and {r 3 ,r 6 }. According to convention, 23 
the a component of the doubly degenerate coordinates is defined such that it is symmetric 
with respect to a <r v reflection plane, whereas the b component is antisymmetric. The 
x, y and z components of the triply degenerate coordinate are also defined according to 

established convention. 24 

The precision of the central difference numerical procedures used to obtain the force 
constants has been closely monitored: the uncertainty in the harmonic frequencies should 
be less than 0.1 cm -1 and in the fundamental frequencies less than 0.5 cm '. The anhar- 
monic analyses have been performed with the SPECTRO program package 25 which uses 
second-order perturbation theory; Ca 3 has been treated as a symmetric top 23 and Car has 
been treated as a spherical top. 25 The coupled-cluster calculations were performed with the 
TITAN set of programs 22 interfaced to the MOLECULE-SWEDEN suite of programs. 25 
The MRCI calculations were performed with MOLECULE-SWEDEN. 

Results and Discussion 

A. Equilibrium structures and binding energies 

Table 1 contains the equilibrium bond distances, rotational constants and dissociation 
energies (atomization energies) of Ca 3 and Ca 4 computed in this study. For Ca 3 the 
CCSD level of theory substantially underestimates the D e value (i.e., 5.7 kcal/mol or 47% 
compared to MRCI) and therefore yields a bond length 0.22 a 0 too long. However, the 
CCSD(T) level of theory exhibits a significant improvement over the CCSD results: D e 
is only 0.9 kcal/mol less than the MRCI result and the equilibrium bond distance differs 
from the MRCI value by only 0.01 o 0 . A similar situation is found with Ca 4 — the CCSD 
level of theory substantially underestimates the binding energy and gives an equilibrium 
bond distance 0.11 a 0 too long. However, as was found for the Be and Mg clusters, 11 ' 12 
the CCSD level of theory performs better for Ca 4 than it does for Ca 3 . The CCSD(T) 
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results for Ca 4 (using the smaller ANO basis set) are in very good agreement with the 
respective MRCI quantities. The CCSD(T) D e is only 0.2 kcal/mol less than the MRCI 
value and the CCSD(T) r e is actually 0.02 a 0 shorter than the MRCI r e . As with our 
earlier studies 11 - 12 of small Be and Mg clusters, it thus appears that the CCSD(T) level of 
theory provides a very good description of electron correlation effects in small Ca clusters. 

It is interesting to note that comparison of the equilibrium bond distances obtained 
with the CCSD, CCSD(T) and MRCI methods suggests that the bonding m small Ca 
clusters is intermediate between the bonding in small Be and small Mg clusters, as would 
be expected based on the bulk binding energies. 29 For Be 3 , where sp hybridization is known 
to play an important role in the binding, the CCSD level of theory yields a reasonable 
equilibrium bond distance when compared to MRCI or CCSD(T). However, for Mgs, where 
the binding is more dominated by dispersion, the CCSD level of theory gives a bond 
distance that is significantly too long (0.55 a 0 ). The discrepancy in the CCSD bond length 
for Ca 3 is between that found for Be 3 and Mg 3 , but is much closer to the value obtained 
for Be 3 than that found for Mg 3 . It therefore seems that sp hybridization is an important 
component of the bonding in Ca 3 , though not as important as it is in the bonding of Be 3 . 

For Ca 4 , only the CCSD and CCSD(T) levels of theory could be used in conjunction 
with the larger ANO basis set since the MRCI procedure would have been prohibitively 
expensive. Using the [6s 5p 2d If] ANO basis set, the CCSD(T) equilibrium bond distance 
for Ca 4 is 0.29 a 0 shorter than the analogous Ca 3 value indicating the increased importance 
of sp hybridization in the bonding as the cluster size becomes larger. The best D e values for 
Ca 3 and Ca 4 obtained in this work are 12.1 kcal/mol and 31.5 kcal/mol, respectively. As 
expected, our best computed D e for Ca 4 is substantially larger than the previously best 
computed value. 14 Based on the fact that the separated-atom limit is described better 
than the molecule, the D e values for Ca 3 and Ca 4 are probably underestimated somewhat. 
We say “probably” because it is not certain whether the effects of core-correlation will 
increase or decrease the D e values. However, based on a recent study of Ca 2 by Dyall and 
McLean, 30 it is likely that the effects of core-correlation will not affect the D e values by 

more than 1-2 kcal/mol. 


B. Vibrational frequencies of Ca 3 

The quadratic force constants and harmonic frequencies of Ca 3 obtained at the 
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CCSD, CCSD(T) and MRCI levels of theory are presented in Table 2. The [6s 5p 2d 1/] 
ANO basis set was used with all of these methods. In view of the underestimation of the 
bond strength at the CCSD level, it is not surprising that the CCSD quadratic force con- 
stants and harmonic frequencies are noticeably smaller than the analogous MRCI values. 
Conversely, the CCSD(T) quadratic force constants and harmonic frequencies are in excel- 
lent agreement with the MRCI quantities, being only 1 cm -1 and 2 cm 1 smaller for the 
a\ and e' modes, respectively. The best computed harmonic frequencies of Ca 3 obtained 
in this work are 95 cm -1 and 85 cm -1 for the a\ and e' modes, respectively. These values 
should be reasonably reliable and are probably low rather than high, assuming that the 
binding energy is still somewhat underestimated. The harmonic frequencies of Ca 3 are 
small not only because the bonding in Ca 3 is relatively weak, but also because of the large 
mass of the Ca atom and the 1 j yjm dependence of the harmonic frequency, where m is 
the reduced mass of the system. 

We have computed cubic and quartic force constants using the CCSD(T) method, 
since the CCSD(T) and MRCI harmonic frequencies are in excellent agreement and the 
former approach is significantly cheaper. The complete set of cubic and quartic force con- 
stants and the resulting anharmonic constants are presented in Table 3. Table 4 contains 
the fundamental vibrational frequencies of Ca 3 determined via second-order perturbation 
theory. In addition, Table 4 presents the infrared (IR) intensity of the e' vibration deter- 
mined using the double harmonic approximation. 

The absolute anharmonic contribution to the vibrational frequencies is small, only 
2 cm -1 for both vibrations, and the percentage effect is about half that observed previously 
for the Mg clusters. Combining the MRCI harmonic frequencies with the anharmonic 
correction obtained with the CCSD(T) method gives 93 cm" 1 and 83 cm" 1 as our best 
estimates for the fundamental vibrational frequencies of the a[ and e' modes, respectively. 
Because the IR intensity of the e’ mode is so small, the best prospect for determining the 
vibrational frequencies experimentally is likely to be an indirect technique such as negative 
ion photoelectron spectroscopy. 


C. Vibrational frequencies of Ca4 

The quadratic force constants and harmonic frequencies of Ca 4 , determined at the 
CCSD, CCSD(T) and MRCI levels of theory, are presented in Table 5. A noteworthy point 
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is that the CCSD harmonic frequencies are in better agreement with the MRCI values than 
was found for Ca 3 . This is especially true for a* and w, where the differences are only 
4 cm" 1 and 8 cm- 1 , respectively. This observation supports our earlier analysis regarding 
the increased importance of 5 p hybridisation, and consequently covalent bonding, in Car 
relative to Ca 3 . The CCSD(T) quadratic force constants and harmonic frequencies are in 
excellent agreement with the respective MRCI quantities. Indeed, <1, and u> 2 only differ 
by 1 cm" 1 and u, differs by less than this. These comparisons suggests that the CCSD(T) 
level of theory is closely approaching the n-particle limit for Car. 

The CCSD and CCSD(T) harmonic frequencies obtained with the larger [6s 5 P 2<f 1/] 
ANO basis set demonstrate the importance of using large one-particle basis sets in order 
to obtain highly accurate harmonic frequencies. The difference between the CCSD(T) 
harmonic frequencies in the two basis sets used is larger than the differences due to cor- 
relation treatment among the small basis results. As with Ca 3 , it is expected that the 
CCSD(T)/[6s 5p 2d If] quadratic force constants and harmonic frequencies should be 

very reliable. 

Table 6 contains the complete cubic and quartic force field of Ca 4 obtained at the 
CCSD(T) level of theory with the larger ANO basis set. The resulting anharmomc con 
stants are given in Table 7, while the fundamental vibrational frequencies and IR intensities 
of Ca 4 are presented in Table 8. As with Ca 3 , the absolute anharmonic corrections for 
the vibrational modes of Ca 4 are relatively small at only 3 cm \ 1 cm 1 and 1 cm for 
u, l5 u >2 and W3 , respectively. The percentage effect on the fundamental frequencies is also 
similar to that found for Ca 3 . The IR intensity of the t 2 vibration is 2.1 km/mol, which 
is substantially larger than the IR intensity of the e' vibration in Ca 3 , but again the best 
method to obtain fundamental frequencies from experiment may be an indirect approach. 


D. Summary of CCSD(T) Results for Small Be, Mg and Ca Clusters 

A summary of the equilibrium structures, vibrational frequencies, infrared intensities 
and binding energies for the alkaline-earth trimers is presented in Table 9 and for the 
tetramers in Table 10. The beryllium and magnesium cluster results are taken from our 
previous studies. 11 ’ 12 In all cases, the results are those obtained with the largest ANO 
basis set used in the particular investigation. For the trimers, Table 9 contains the MRCI 
equilibrium bond distance, harmonic frequencies and dissociation energy. The fundamental 
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frequencies were obtained by adding the CCSD(T) anharmonicity to the MRCI harmonic 
frequencies. The IR intensities were determined with the CCSD(T) method. For the 
tetramers, all of the results were obtained at the CCSD(T) level of theory since it was not 
possible to use the MRCI method in conjunction with the larger ANO basis sets for the 
tetramers. Thus the values summarised in Tables 9 and 10 represent the best computed 
quantities to date for the alkaline-earth trimers and tetramers. 

Examination of the binding energies in Table 9 indicates the expected trend based on 
bulk binding energies - that is the binding energies decrease in the order Be, > Ca, > Mg,. 
The vibrational frequencies, on the other hand, decrease in the order Be, > Mg, > Ca,, 
but this is determined in large part by mass effects, since the symmetry internal coordinate 
quadratic force constants for Mg, are smaller 12 than those for Ca,. The IR intensity of 
the e ; mode is small for all of the alkaline-earth trimers. 

Comparison of the D. values of the tetramers shows the same trend as observed 
for the trimers, although the ratios are somewhat different. The binding energy of Be 4 
is significantly larger than that of Be„ though the Be, equilibrium bond distance is only 
0.28 a, shorter than the Be, value. The ratio D t (M,)/B.(M,) is largest for M = Mg 
and it is therefore not surprising that Mg, exhibits the largest reduction in bond length 
(0.50 oo) relative to the trimer. The reduction in the Ca, bond length (relative to the 
trimer) is about the same as that observed for Be even though the ratio Z>«(M,)/.D,(M,) 
is significantly larger for M = Be than for M = Ca. This last observation is probably due 
to the fact that the valence orbitals of Ca are larger than those of Be. 

The harmonic frequencies of the tetramers exhibit the same trend as observed for the 
trimers. Again, the Mg, vibrational frequencies are higher than the analogous Ca, values 
because of the mass effect. The fundamental frequency n for Be, has a large positive 
anharmonic correction that is not observed in either Mg, nor Ca,. This Be, phenomenon 
was explained in some detail previously, 11 and relies on a symmetry argument applicable 
to tetrahedral X, species. Its apparent inapplicability to Mg, and Ca, is probably due 
to several factors, including the much weaker bonds present in Mg, and Ca, relative to 
Be, The IR intensities of the tetramer 1, vibrations are much larger than those calculated 
for the trimers. Nevertheless, the Mg, and Ca, intensities remain very small and it is 
likely that the best prospect for experimental determination of these frequencies will be an 
indirect method like negative ion photodetachment. On the other hand, the IR intensity 
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of the h mode of Be 4 is certainly large enough to allow direct experimental observation 
provided that an experiment can be designed which will produce enough Be 4 . 

Conclusions 

The CCSD, CCSD(T) and MRCI electron correlation methods have been used to 
investigate the equilibrium structures, vibrational frequencies and binding energies of the 
Ca 3 and Ca 4 metallic clusters. In agreement with our earlier studies of the analogous Be 
and Mg clusters, it is found that the CCSD(T) method reproduces the MRCI r«, harmonic 
frequencies and D e values very well. The agreement between these methods is somewhat 
better for Ca 4 , but is still very good for Ca 3 . Complete cubic and quartic force fields 
of Ca 3 and Ca 4 have been determined with the CCSD(T) method m conjunction with a 
large ANO basis set and have been used to evaluate the anharmonic corrections needed to 
compute the fundamental frequencies. The anharmonic corrections have been determined 
via second-order perturbation theory. The absolute value of the anharmonic corrections 
is relatively small, although as a percentage relative to the fundamental frequencies they 
are similar to those observed previously for the Mg clusters. In spite of the fact that Ca is 
larger and more polarizable than Mg, and that Ca clusters are more strongly bound than 
Mg clusters, the IR intensities of the e' mode of Ca 3 and of the t 2 mode of Ca 4 are small 
and similar to the analogous Mg quantities. It is unlikely that direct observation of these 

fundamentals will be possible. 

The MRCI and CCSD(T) equilibrium structures, vibrational frequencies and binding 
energies of the alkali metal (Be, Mg and Ca) trimers and tetramers have been summarized. 
The binding energies of the trimers and tetramers follow the bulk metal binding energies 
although the ratios of the small cluster D e values do not agree with the bulk metal ratios. 
The vibrational frequencies follow a different trend as the Mg n (n=3,4) frequencies are 
larger than the respective Ca values, but this is due to the larger mass of the Ca atom, 
since the symmetry internal coordinate force constants (which are independent of mass) 
for the Ca trimer and tetramer are larger than the respective Mg quantities. 
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Table 1 


Total energies (Eh), bond lengths (a 0 ), rotational 
constants (MHz) and binding energies (kcal/mol) for Ca 3 and Ca 4 . 




Basis 

E a 

r e 

A 

B 

D e 

D o 

Ca 3 

CCSD 

[6s 5 p 2d If] 

0.370728 

8.097 

1378 

689 

6.4 

6.1 


CCSD(T) 

[6s 5 p 2d If] 

0.378316 

7.884 6 

1453 

726 

11.2 

10.8 


MRCI 

[6s 5 p 2d If] 

0.379234 

7.874 

1457 

728 

12.1 

11.7 

C^4 

CCSD 

[5s 4p Id] 

0.168123 

7.935 

717 


14.3 

13.6 


CCSD(T) 

[5s 4p Id] 

0.180604 

7.806 

741 


22.2 

21.4 


MRCI 

[5s 4 p Id] 

0.179059 

7.824 

738 


22.4 

21.7 


CCSD 

[6s 5 p 2d If] 

0.180579 

7.694 

763 


20.9 

20.1 


CCSD(T) 

[6s 5 p 2d 1/] 

0.197496 

7.591 c 

784 


31.5 

30.6 


a The energy for Ca 3 is reported as —(£+2030) and for Ca 4 as —(£+2707). 
b Vibrationally averaged bond lengths: r g = 7.917 ao and r a = 7.915 ao- 
c Vibrationally averaged bond lengths: r g = 7.615 ao and r a = 7.613 ao- 
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Table 2 

Symmetry internal coordinate force constants (aJ/A ) 
and harmonic frequencies (cm 1 ) for Ca 3 . 



F n 

F 22 



CCSD 

0.04415 

0.08141 

75 

72 

CCSD(T) 

0.06952 

0.10835 

94 

83 

MRCI 

0.07126 

0.11460 

95 

85 
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Table 3 

Non-zero cubic (aJ/A 3 ) and quartic (aJ/A 4 ) force 
constants and the anharmonic constants (cm -1 ) for Ca 3 . 


■Fin 


-0.133 

F\2a2a = 

F 12626 

-0.152 

F 2a2a2a — 

— F 2a 2b2b 

-0.105 

-Pllll 


0.156 

F\\2a2a = 

Fli2b2b 

0.163 

F 12a2a2a = 

~Fj2a2b2b 

0.116 

F2a2a2a2a ~ 

= F 2 b2b2b2b = 3-P 2a2a2b2b 

0.239 


Anharmonic Constants 

*ii 

-0.60 

X 2 1 

-1.38 

%22 

-0.41 

9 22 

0.18 
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Table 4 

Comparison of the CCSD(T) harmonic and fundamental frequencies 
of Ca 3 (cm -1 ). Infrared intensities (km/mol) are also included. 


Mode 

U> 

V 

u - V 

I 


94 

92 

2 

0 

e' 

83 

81 

2 

0.4 
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Table 5 

Symmetry internal coordinate quadratic force constants (aJ/A 2 ) 
and harmonic frequencies (cm -1 ) for Ca 4 . 



Basis 

F n 

F 22 

F 33 

^i(ai) 

u> 2 (e) 

<^ 3 (^ 2 ) 

CCSD 

[5s 4 p Id] 

0.05350 

0.12546 

0.08688 

95 

73 

86 

CCSD(T) 

[ 5 s 4 p Id] 

0.07166 

0.14253 

0.10439 

110 

78 

94 

MRCI 

[ 5 s 4 p Id] 

0.06960 

0.14091 

0.10397 

109 

77 

94 

CCSD 

[ 6 s 5 p 2 d 1 /] 

0.07463 

0.15725 

0.11282 

113 

82 

98 

CCSD(T) 

[ 6 s 5 p 2 d 1 /] 

0.09446 

0.17250 

0.13062 

127 

86 

105 
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Table 6 

Non-zero cubic (aJ/A 3 ) and quartic (aJ/A 4 ) force 
constants for Ca4. 


-Pill 


-0.110 

F 12a2a = ^12626 


-0.140 

F\ZxZx = FiZyZy = FiZzZz 


-0.126 

-p2a2a2a = — -p2a2626 


-0.074 

F2aZzZz ~ — 2i^2a3z3z = “ 

■2F2a3y3y = ^-^263x33: — ^-F263y3y 

-0.190 

^3a;33/3z 


-0.013 

-Fun 


0.077 

Pll2a2a = -Pli2&2& 


0.091 

■Pll3x3i — -Fil3y3y — -^1 1 3 z 3 z 

0.076 

F\2a2a2a — “‘Fl2a2&2& 


0.058 

Pl2a3z3z = 2 -F\ 2 a 3 x 3 x = 

~2Fi2a3y3y = ~^Fi2b3x3x — ^|‘-?12b3y3y 

0.119 

-P\ 3x3y3z 


-0.001 

jp2 a 2 a 2 a 2 a = -?2b262&26 — 3i ? 2a2a2626 

0.162 

-p2a2a3z3z 


0.176 

F2b2b3z3z 


0.044 

F2a2a3y 3y — -F2a2a3x3x 

z j(F2a2o3z3z + 3-p2*263z3z) 

0.077 

a F26263y3t/ = p26263z3a; = 

^(3i 7 2a2a3z3z + F26263z3z) 

0.143 

F2a263y3y — -F2a263x3x 

= ^■(P2a2a3z3z “* p2&263z3z) 

0.057 

F^xZxZxZx — F^yZyZyZy 

FzzZzZzZz 

0.262 

F3x3x3y3y — F^xZxZzZz 

FzyZyZzZz 

0.043 


a Dependent force constants related to F^a2aZzZz and -p 262 !> 3 z 3 z- 
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Table 7 

Anharmonic constants (cm -1 ) for Ca 4 . 


Anharmonic Constants 


*11 

-0.31 

*21 

-0.52 

*22 

-0.13 

®31 

-0.83 

^32 

-0.35 

£33 

-0.18 

922 

0.09 

9 33 

0.03 

^23 

-0.05 

^33 

-0.03 
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Table 8 

Comparison of the CCSD(T) harmonic and fundamental frequencies 
of Ca 4 (cm -1 ). Infrared intensities (km/mol) are also included. 


Mode 

U) 

V 

U) — V 

I 


127 

124 

3 

0 

e 

86 

85 

1 

0 

^2 

105 

104 

1 

2.1 
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Table 9 

Summary of results for the alkaline earth tnmers. 


Be 3 4.200 
Mg 3 6.373 
Ca 3 7.874 


w i( a :i) u, 2(e') v \ 


427(0.5) 469 410 22.5 

115(0.2) 101 109 6.3 

85(0.4) 93 83 12.1 


• units are »„ for r„ cm' 1 for the harmonic and fundamental frequencies, and kcal/mol 
for D. The value in parentheses is the IR intensity in km/mol. The Be, results are from 
reference 11 and the Mg, results are from reference 12. See text for details of correlate 

methods used. 
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Table 10 

Su mm ary of CCSD(T) results for the alkaline earth tetramers. a 



r e 

wi(ai) 

a> 2 (e) 

<^3(^2) 


^2 

^3 

D e 

Be4 

3.921 

663 

469 

571(29.7) 

639 

455 

682 

79.5 

Mg 4 

5.877 

192 

147 

171(2.4) 

184 

143 

167 

23.9 

Ca 4 

7.591 

127 

86 

105(2.1) 

124 

85 

104 

31.5 


° Units are a 0 for r e , cm -1 for the harmonic and fundamental frequencies, and kcal/mol 
for D e . The value in parentheses is the IR intensity in km/mol. The Be 4 results are from 
reference 11 and the Mg 4 results are from reference 12. 
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